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Abstract. The magnetic behavior of thin ferromagnetic itinerant-electron films is 
investigated within the strongly correlated single-band Hubbard model. For its 
approximate solution we apply a generalization of the modified alloy analogy (MAA) 
to deal with the modifications due to the reduced translational symmetry. The 
theory is based on exact results in the limit of strong Coulomb interaction which 
are important for a reliable description of ferromagnetism. Within the MAA the 
actual type of the alloy analogy is determined selfconsistently. The MAA allows, 
in particular, the investigation of quasiparticle lifetime effects in the paramagnetic 
as well as the ferromagnetic phase. For thin fcc(100) and fcc(lll) films the layer 
magnetizations are discussed as a function of temperature as well as film thickness. The 
magnetization at the surface-layer is found to be reduced compared to the inner layers. 
This reduction is stronger in fcc(100) than in fcc(lll) films. The magnetic behavior 
can be microscopically understood by means of the layer-dependent spectral density 
and the quasiparticle density of states. The quasiparticle lifetime that corresponds to 
the width of the quasiparticle peaks in the spectral density is found to be strongly 
spin- and temperature-dependent. 
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1. Introduction 



Remarkable advances in thin film technology have recently led to active interest in 
the nature of magnetism in ultrathin films, at surfaces and multilayer structures. The 
influence of the reduced dimensionality on the magnetic behavior of 3d transition metals 
has been extensively studied both experimentally |], §, ||, |[ ||, ||, [7] and theoretically 
fl, |, [II], [12], [H| [14|, [L5|, [U|. On the experimental side it was shown that ultrathin 
transition metal films can display long range ferromagnetic order from a monolayer on 



g. Although the Mermin- Wagner theorem [17| requires the transition temperature 
to vanish for perfectly isotropic two dimensional systems, it was shown theoretically 
that even a small amount of anisotropy may lead to magnetic order with a substantial 
transition temperature |18|, |19|, ED . In real materials magnetic anisotropy is always 



present by virtue of either the dipole interaction or the spin-orbit coupling. 

Theoretically, the T = K properties of thin transition metal films have been 
addressed by ab initio calculations within the density functional theory in the local 
density approximation & E| [ 



11 



However, these approaches are strictly based on 
a Stoner-type model of ferromagnetism and, therefore, treat electron correlation effects 
which are responsible for the spontaneous magnetic order on a low level. In addition 
they are restricted to ground state properties only. To overcome this restriction, for 
example, a generalization of the fluctuating local moment method has been used 
to calculate the temperature-dependent electronic structure of thin ferromagnetic films. 
However, the layer magnetizations at finite temperatures and the magnetic short range 
order are needed as an input. In Ref. |T2[ magnetic phase transitions in thin films 



are investigated via a mapping of the ab initio results onto an effective Ising model. 
Hasegawa calculates the finite temperature properties of thin Cu/Ni/Cu sandwiches by 
use of the single-site spin-fluctuation theory . 



For the understanding of the thermodynamical properties of thin film magnetism 
theoretical investigations on rather idealized model systems have proven to be a good 
starting point. In this context several authors have focused on localized spin models 
like the Heisenberg model [22, 23, 24, 25 1. For example, the mechanism that leads 



to the experimentally observed temperature induced reorientation of the direction of 
magnetization in thin Fe and Ni films || |7| was investigated in great detail [^, |27|j . 
On the other hand it is by no means clear to what extent the results obtained by 
localized spin models are applicable to transition metal films, where the magnetically 
active electrons are itinerant. 

The aim of the present paper is to study the interplay between strong electron 
correlations and the reduced translational symmetry due to the film geometry within 
an itinerant electron model system. In particular we are interested in the influence of 
the reduced dimensionality on spontaneous ferromagnetism and the spin-, layer- and 
temperature-dependent electronic structure. For this purpose we restrict ourselves, at 



present, to the investigation of the single-band Hubbard model ||28|| , which includes the 
minimum set of terms necessary for the description of itinerant-electron magnetism. 
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The Hubbard model was originally introduced to explain band magnetism in transition 
metals and has become a standard model to study the essential physics of strongly 
correlated electron systems over the years. It is clear that a realistic and quantitative 
description of ferromagnetism in transition metals requires the inclusion of the 
degeneracy of the 3d-bands [^, ^0], |3T], [32], [33], |3|, |35|. Although the band-degeneracy is 
neglected in our model study, we believe that a treatment of electron correlation effects 
well beyond Hartree-Fock theory will provide important insight into generic properties 
of thin film ferromagnetism. For example, contrary to the expectation on the basis of the 
well-known Stoner criterion, the magnetic order at the film surface may be reduced and 
less stable compared to the inner layers if electron correlations are taken into account 
properly f36f . 

Despite its apparent simplicity no general solution exists until now for the Hubbard 
model. However, recently exact results have been obtained by finite temperature 
quantum Monte Carlo calculations in the limit of infinite dimensions |3~7 , |3~0f which prove 
the existence of ferromagnetic solutions for intermediate to strong Coulomb interaction 
U. In addition the decisive importance of the lattice geometry, i.e. the dispersion and 
distribution of spectral weight in the non-interacting (Bloch) density of states (BDOS), 
on the magnetic stability was stressed by several authors J3(| [37], [38], [E], |41], £42|. A 
reasonable treatment of electron correlation effects led to an argument for the stability 
of ferromagnetism which is decisively more restrictive the well known Stoner criterion. 
A BDOS with large spectral weight near one of the band edges is an essential ingredient 
for ferromagnetism. The thermal stability of ferromagnetic solutions is favored by a 
strong asymmetry in the BDOS This behavior of the BDOS is found, for 

example, in non-bipartite lattices like the fee lattice. 

Due to the broken translational symmetry even more complications are introduced 
to the highly non-trivial many-body problem of the Hubbard model. Thus we require an 
approximation scheme which is simple enough to allow for an extended study of magnetic 
phase transitions and electronic correlations in thin films. On the other hand it should 



be clearly beyond Hartree-Fock (Stoner) theory which has been applied previously [14 



since we believe a reasonable treatment of electron correlation effects to be vital for a 
proper description of ferromagnetism especially for non-zero temperatures. 

In this context interpolating theories which are essentially based on exact results 
obtained by the 1/U perturbation theory first introduced by Harris and Lange |44| ^ 



have proven to be a good starting point p3| . A theory that reproduces the rigorous 
strong coupling results in a conceptual clear and straightforward manner is given by the 
spectral density approach (SDA) which has been discussed with respect to spontaneous 
magnetic order for various three-dimensional [H6|, [I7|, |4^| as well as infinite dimensional 
fll?| , lattices. A similar approach applied to a multiband Hubbard model led to 
surprisingly accurate results for the magnetic key quantities of the prototype band 



ferromagnets Fe, Co, and Ni |^3[. A generalization of the SDA to systems with reduced 
translational symmetry has recently been given in Refs. |T5|, |T^], which led, for example, 
to the description of the temperature-driven reorientation transition within an itinerant- 
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electron film |16fl . However, a severe limitation of the SDA results from the fact that 
quasiparticle damping is neglected completely. To tackle this problem a modified alloy 
analogy has been proposed |48], |4£j which is closely related to the SDA but includes 



quasiparticle damping effects in a natural way. For bulk systems it was found that 
the magnetic region in the phase diagram is significantly reduced by the inclusion of 
damping effects. By comparison [^, [5U| with exact results for the fee lattice in the limit 
of infinite dimension and intermediate Coulomb interaction it is clear that the Curie 
temperatures are somewhat overestimated within the MAA. However, the qualitative 
behavior of the ferromagnetic solutions and in particular the dependence of the Curie 
temperature on the band occupation is found to be in good agreement with the exact 
results. 

In the present work we want to apply the MAA to systems with reduced 
translational symmetry. For this purpose the paper is organized in the following way: In 
the next section we will give a short introduction to the underlying many-body problem. 
The concept of an alloy analogy for the Hubbard film is developed in Sect. [| In Sect. f| 
we will generalize the MAA to systems with reduced translational symmetry. The results 
of the numerical evaluations will be discussed in Sect. § in terms of temperature- and 
layer-dependent magnetizations, the quasiparticle bandstructure and the quasiparticle 
densities of states. We will end with a conclusion in Sect. ||. 

2. The many-body problem of the Hubbard film 

Let us first introduce the notation used to deal with the film geometry. Each lattice 
vector of the film system is decomposed into two parts according to: 

R<ia = R-i + r a . (1) 

Rj denotes a lattice vector of the underlying two-dimensional Bravais lattice with N 
sites. To each of theses lattice sites there is associated a rf-atom basis r a (a — 1, . . . , d) 
which refers to the d layers of the film. The same labeling with Latin and Greek indices 
applies for all quantities related to the film geometry. Within each layer we assume 
translational invariance. Then a two-dimensional Fourier transformation with respect 
to the Bravais lattice can be applied. 

Using this notation the Hamiltonian for the single-band Hubbard film reads: 

W = ( T tf ~ ^ijViaaCjP* + n iacT n ia - a . (2) 

i,j,a,[3,cr i,a,a 

Here Q acr (cJ ao .) stands for the annihilation (creation) operator of an electron with spin 
a at the lattice site Rj a , n iaa = c\ acr c iaa is the number operator. U denotes the on-site 
Coulomb matrix element and /i the chemical potential. T?j is the hopping integral 
between the lattice sites Rj a and Rj/3. A two-dimensional Fourier transformation yields 
the corresponding dispersions 

T f = -J2 T^e-^-^\ (3) 
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Here and in the following k denotes a wave- vector from the underlying two-dimensional 
(surface) Brillouin zone. Further we define T 0a = T? a = -k J2kT£ a — const, which gives 
the center of gravity of the a-th layer in the BDOS. 

The basic quantity to be calculated is the retarded single-electron Green function 

G°f a (E) = ((c iaa ;c] 0a )) E . (4) 

From Gfj a (E) we can obtain all relevant information on the system. After a two- 
dimensional Fourier transformation one obtains from Gff a (E) the spectral density 

S£{E) = -hmG%(E), (5) 

7T 

which represents the bare lineshape of a (direct, inverse) photoemission experiment. 
The diagonal elements of the Green function determine the spin- and layer- dependent 
quasiparticle density of states (QDOS): 

PUB) = ^ E S£(E ~ A*) = -\^GZ{E - n). (6) 

k 

Via an energy integration one immediately gets from p aa (E) the band occupations 

oo 

n aa = (n iaa ) = J dEf^(E)p aa (E). (7) 

— oo 

(. . .) denotes the grand-canonical average and f-{E) is the Fermi function. Here the site 
index % has been omitted due to the assumed translational invariance within the layers. 
Ferromagnetism is indicated by a spin-asymmetry in the band occupations n aa leading 
to non-zero layer magnetizations m a = n a ^ — n a i- The mean band occupation n and 
the mean magnetization m are given by n = \ J2 acr n aa and m = \ J2 a m «i respectively. 
The equation of motion for the single-electron Green function reads: 

E P + - IT - Kl(E)} G]f a (E) = hStf. (8) 

h 

Here we have introduced the electronic self-energy T,^- (T (E) which incorporates all effects 
of electron correlations. 

For later use we want to define the moments of the Green function 



oo 

— Im J dEE m G?f a (E). (9) 



The usefulness of the moments M^ a/S (m = 0,1,2,...) results from the fact that 
an alternative but equivalent representation can be derived by use of the Heisenberg 
representation of the creation and annihilation operators. Thus M^ al3 can be 



calculated up to the desired order m directly from the Hamiltonian (|2|) itself |46|, |5l| : 

= ([ [...[c ta(T ,ni...,ni , 4J + ). (io) 

m-times 

Here [...,.. •]-(+) denotes the commutator (anticommutator). Eqs. @ and (0) impose 
rigorous sum rules on the Green function and the self-energy which have been recognized 
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to state important guidelines when constructing approximate solutions for the Hubbard 
model P3). For example, the high energy expansion of the Green function is directly 



determined by the moments M^ a ^ . It has been shown El that the sum rules are 



especially important in the limit of strong Coulomb interaction: Being consistent with 
the sum rules up to the order m = 3 states a necessary condition in order to reproduce 
the exact results of the l/£7-perturbation theory Furthermore, the m = 3 sum 



rule turns out to be of particular importance what concerns the stability of ferromagnetic 



solutions in the Hubbard model 43 



The sum rules up to order m = 3 will be exploited in Sect. f| for the construction 
of a modified alloy analogy (MAA) to the Hubbard film. First we want to introduce the 
concept of the alloy analogy approach for systems with reduced translational symmetry. 

3. The alloy analogy concept for the Hubbard film 



The main idea of the conventional alloy analogy approach |52] is to consider, for the 



moment, the — cx-electrons to be "frozen" and to be randomly distributed over the sites 
of the lattice. Then a propagating a-electron encounters a situation which is equivalent 
to a fictitious alloy: At empty lattice sites it finds the atomic energy Ei a , at sites 
with a — a-electron present the atomic energy E 2a . These energy levels are randomly 
distributed over the lattice with concentrations x\ a and x 2a which correspond to the 
probabilities for the cx-electron to meet these respective situations. Note that at this 
point it is not at all clear what choice of the energy levels and concentrations gives the 
best approximation for the initial Hamiltonian. However, an "optimal" choice of the 
alloy analogy parameters should by some means account for the itineracy of the — o- 
electrons (see Sect. §). In the present film system the energy levels and concentrations 
may, in addition, exhibit a layer- dependence. Thus the alloy analogy for the Hubbard 
film is described by 4 • d a priori unknown parameters 

E&\ x£\ E%\ xg a = l,...,d (11) 

For the solution of the fictitious alloy problem given by ([0]) the coherent potential 
approximation (CPA) |53) provides a well known method. The CPA has been realized to 



be the rigorous solution of the alloy problem in the limit of infinite dimensions |54j where 
the single-site aspect used in the derivation of the CPA becomes exact. In this sense the 
CPA can be termed to be the best single-site approximation to the alloy problem. Due 
to the single-site aspect and the assumed translational invariance within the layers we 
have Y%P{E) = 5fPY, aa {E). The implicit CPA equation |53j for the self-energy is readily 



formulated via an effective medium approach similar to the one discussed in Ref. [|5Tf : 

0= T x ( Z ] ~ Eaa ^ ~ T ° a (12) 

p=i,2 1 — \Gf£[Ep$ — E aa (E) — T 0a ] 

In addition the self-energy appears implicitly in the expression for the local Green 
function which is given by matrix inversion from (Rt) after applying a two-dimensional 
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Fourier transformation: 



h 



( E+n-T^-Vx^E) 



rp\2 



iV k 



V 



Eqs. ( p"2D and (fl3|) have to be solved selfconsistently to obtain Gf£(E) and Ti aa (E). 



(13) 



a 3 



4. The modified alloy analogy 

Up to now nothing has been said about the actual choice of the atomic energies E$ and 
the corresponding concentrations (p = 1, 2). In the conventional alloy analogy (AA) 



|52| the alloy parameters (|IT|) are taken from the zero-bandwidth limit which directly 
corresponds to the assumption of strictly "frozen" — cr-electrons: 



E\ a 

~(«) 
x la ■ 



1 n a — a 



5(a) 
2g 

~(a) 
X 2a ■ 



T 0a + U 



n r 



(14) 



However, it was soon realized that the AA is not able to describe itinerant 
ferromagnetism |55|| . This is closely related to the fact that the energy levels E$ arc 
rigid and, in particular, spin-independent quantities within the AA. Further it is known 
H| §§ that the AA fulfills the sum rules ©, (0) up to the order m = 2 only and fails 
to reproduce the correct strong coupling behavior. Note that within the AA the energy 
levels |W are layer-independent (for uniform To a ) which is a crude approximation 
since a possible layer-dependence in the quasiparticle spectrum is suppressed almost 
completely. 

The basic idea of the MAA is to exploit the information provided by the non- 
trivial but exact results in the limit of strong Coulomb interaction (U/t ^> 1) H3, [IIJ 
to determine the energy levels and concentrations (|TT|). This can most elegantly be 
achieved by imposing the sum rules (||), (|T0|) on the CPA equation (|12"D |50| , [43| : By 
inserting the high energy expansion of the self-energy Yi aa (E) and the local Green 



function G°**(E), which are determined by the sum rules, the CPA equation ([12]) can be 
expanded in powers of 1/E. Taking into account the sum rules up to the order m = 3 
unambiguously determines the parameters . Then the exact strong coupling 

results are reproduced automatically j£| • Note that due to the single site aspect of the 
CPA only the local terms of the 1/U perturbation theory are reproduced. On the other 
hand the MAA is not restricted solely to the strong coupling limit but is also applicable 
for intermediate interaction strengths where it has an interpolating character j|8], |49|j . 
Following this procedure yields the energy levels and concentrations of the MAA for the 
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Hubbard film: 



E {a) 

pa 



T 0a + U + B a J±) 



x 



(a) 

1(7 



(U + B a _ a - T 0a ) 2 + 4Un a _ a (T 0a - B a _ a ) 
B Q _ a + U(l - n a _ a ) - E$ 



^2a ~ ^la 



(15) 



T («) 

X 2cr 



1 - T (a) 

1 x la 



An alternative derivation of the MAA for bulk systems which is based on physical 
arguments can be found in Refs. p8| , fi9| . Note that the expressions for E^> and in 
(|j~5|) are directly related to the position and the weight of the two poles of the spectral 
density within the SDA. Eqs. (O) are obtained from the SDA results ESI, H9[ if the 



electron dispersion is replaced by the center of gravity of the non-interacting band. The 
energy levels and concentrations ([15]) are not only dependent on the model parameters 
T 0a and U but also on the band occupations n a -a and the so-called bandshift B a _ a . The 



bandshift that is introduced via the fourth moment M-f}°f consists of higher correlation 
functions: 



B a - 



1 



■ Oo" 



n, 



E T t ?(cl 



Cj(3- a {2nia<j-l)). 



(16) 



-ct(1 ' ' ' n-a i , j 

Nevertheless -B Q _ CT can exactly be calculated [^, ^| by use of the local Green function 
and the self-energy: 



B a - 



1 



n a -M 



+oo 

l -hnjdEf_{E) ^ a _ a (E-fi) 



— CO 



[{E - E a ^(E -fi)- T 0a )G?«a(E -n)-h] (17) 

In the strict zero-bandwidth limit B a _ a is identical to To Q and the MAA ( |T5D reduces 
to the conventional alloy analogy (P^). However, as soon as the hopping is switched on, 
the bandshift B a _ a , which is for strong Coulomb interaction proportional to the kinetic 
energy of the — cr-electrons in the a-th layer [0, has to be calculated selfconsistently 



by iteration. Thus, via (|15D the type of the underlying alloy changes in each step of the 
iteration process. In this sense B a _ a accounts for the itineracy of the — cx-electrons. In 
the paramagnetic phase there are only minor differences in the quasiparticle spectrum 
between MAA and AA. However, the bandshift may get a real spin-dependence for 
special parameter constellations. Thus B a _ ff may generate and stabilize ferromagnetic 
solutions which are excluded within the AA. It is worth to stress that the energy levels 
and concentrations are implicitly temperature-dependent via n a - a and B a _ a leading, 
therefore, to a temperature dependent electronic structure. 

The evaluation of the MAA requires the solution of two nested selfconsistency 
cycles. One starts with an initial guess for the band occupations n a _ a and the bandshift 
B a _ a which determine the energy levels and concentrations (0). Via the CPA-equation 
(|T2|) and ( |T3D the corresponding self-energy and Green function can be calculated 
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selfconsistently. With this solution new values for n a _ cr and B a _ a are obtained via 
(0) and (|ITD . This procedure is iterated until convergence is achieved. For efficiency 
reasons the numerical evaluations of the integrals in (|7|) and (17) are performed via 



discrete Matsubara sums on the imaginary energy axis [pD[ . Only the spectral density 
and the quasiparticle density of states are calculated on the real axis at the end of each 
selfconsistency procedure. 

5. Results and Discussion 

For the numerical evaluations we consider in the present work thin fee films with an 
(100) as well as an (111) surface and a film thickness up to d = 15. The hopping integral 
between nearest neighbor sites is chosen to be uniform throughout the film and is set 
to t = — 0.25eV. All other hopping integrals as well as To a are set to zero. For an fee 
bulk system this yields a total bandwidth P1/ bulk = 4eV of the non-interacting system. 
Further, we keep the on-site Coulomb interaction fixed at U = 50 eV which clearly refers 
to the strong coupling regime. In all calculations the total band occupation is kept fixed 
at the representative value n = 1.6. Bulk calculations within the MAA have shown 



that for the fee lattice ferromagnetic order is possible for all band occupations above 
half filling n > 1. 



Table 1. Number of nearest neighbors (n.n.) for the fee lattice as well as the fcc(lOO) 
and the fcc(lll) film structure. In addition the moments (see Eq. (19)) of the 

BDOS are given. A« refers to the surface layer, while the moments of all other layers 
(2 < a < d — 1) are equal to the respective bulk values A^K 



bulk 


fee 


n.n 


12 


A?> 


12 


Af 


-48 


film 


(100) 


(111) 


+ 1 


4 


3 


n.n. 


4 


6 


-1 


4 


3 


Ai 2) 


8 


9 




-24 


-30 




0.667 


0.750 


A (3) /A (3) 


0.5 


0.625 



For both film structures considered here, fcc(100) and fcc(lll), all nearest neighbors 
are placed in the same or in the adjacent layer. The number of nearest neighbors in 
these two film geometries are given in Tab. [|. The corresponding dispersions can 
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then be written as: 

!Toa + tj\\(\s), a = P 

t 7 ±(k), a = P-l 

t(T±(k))*, c* = /3 + l. 

According to (H) one gets 



(loo) /■. \ o r f^x + ky k x k y , 
< ]r > - ' cos( -) + cos( 1 



7|| (k) 

7 [ 100) (k) = 1 + t '-r- + f '— - + e" 
for the fee (100) film geometry and 

7 | ( | m) (k) = 2 [(308(^/3/8^ + v/l/8^) + cos(y/3/8k x - + 



cos(yl/2A; y ) 

7 (_ in) (k) = 1 + e-^Vs/^+VVskv) + g-KV 3 / 8 ^-"^ 1 /^) 

for fcc(lll). Here the lattice constant is set to a = 1. The layer-dependent Bloch 
density of states po a (E) = p a a(E)\u=o for a five layer film is plotted in Fig. [I] for both 
film structures considered. The BDOS is strongly asymmetric and shows a distinct layer 
dependence. Considering the moments 

oo 

Ai m) = — / dE(E - T 0a r Poa (E) (19) 

— oo 

of the BDOS yields that the variance A^ 2 ) as well as the skewness are reduced at 
the surface layer compared to the inner layers due to the reduced coordination number 
at the surface (see Tab. [I]). 

The charge distributions n a as well as the layer magnetizations m a are determined 
by the selfconsistently calculated QDOS © via (0). The chemical potential \i and the 
band centers T 0a are assumed to be uniform throughout the film, allowing, therefore, 
for charge transfer between the layers. However, in the actual calculation the difference 
in the occupation numbers n a turns out to be very small (< 3%). 

The layer-dependent magnetizations m a together with the mean magnetization m 
for a five layer film are plotted in Fig. ^| as a function of temperature for both film 
structures. With respect to the overall shape the magnetization curves show the usual 
Brillouin-type behavior. However, the surface magnetization is found to be reduced 
compared to the inner layers for all temperatures. The reduction is particularly strong 
for the fcc(lOO) film geometry and leads to a non-saturated groundstate whereas the 
fcc(lll) film is fully polarized at T = 0. The enhanced surface effects in the fcc(lOO) 
structure are related to the higher percent of missing nearest neighbors at the surface 
layer which is 1/3 for fcc(lOO) and 1/4 for fcc(lll). 

We want to emphasize that the finding of a reduced surface magnetization cannot 
be explained by the well-known Stoner criterion of ferromagnetism. Since the variance 
of the BDOS is reduced at the surface layer (A^ < A^, see Tab. |I|) due to the reduced 
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Figure 1. Layer-dependent Bloch density of states po a (E) of a five layer film for (a) 
the fcc(100) and (b) the fcc(lll) structure, a = 1 denotes the surface layer, a = 3 the 
central layer. Further parameters: t = — 0.25eV, Tq q = 0. 



coordination number one might intuitively expect the magnetization at the surface to be 
more robust than in the bulk. However, as discussed in Sect. [I], intensive investigations 
of strongly correlated electron systems well beyond Hartree-Fock (Stoner) theory clearly 
point out the importance of a large skewness A^ 3 ^ for the stability of ferromagnetism 
p9| , f40| , (43]. Since the skewness of the BDOS is strongly reduced at the surface (see 
Tab. |I|) this explains the trend of a reduced surface magnetization. The above argument 
can be checked by considering the BDOS of the surface layer (see. Fig. |IJ) as an input for 
an additional MAA calculation. Doing so we find that a "separated" surface layer would 
be ferromagnetic for the fcc(lll) but paramagnetic for the fcc(100) film structure. In 
this sense the surface layer of an fcc(100) film is magnetized only because of the effective 
field induced by the ferromagnetically ordered inner layers. 

The Curie temperature is found to be unique for the whole film. Note, that 
although the mean magnetization is reduced for the fee (100) film with respect to 
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Figure 2. Layer magnetizations m a as a function of temperature for (a) an fcc(100) 
and (b) an fcc(lll) five layer film, a = 1 corresponds to the surface layer, a = 3 to 
the central layer. Further parameters: U = 50 eV, t = — 0.25eV, n = 1.6. 



fcc(lll) the corresponding Curie temperature is enhanced (Tq (d = 5) = 1140 K, 
T { c ll \d = 5) = 1050 K). The inner layers that are fully polarized at T = K for both 
film structures appear to be magnetically more stable for fcc(100) compared to fcc(lll). 
Again, this trend can also be seen in an additional MAA calculation for the BDOS of 
the respective central layers. The Curie temperatures converge to the corresponding 
bulk value (T c bulk = 1050 K) for rf( m ) w 3 and w 6. 

In Fig. the surface-, center-, and mean magnetization (m 1; m c and m) are shown 
as a function of the film thickness. The surface magnetization is reduced compared to 
the mean magnetization. The reduction is weak for fcc(lll) films but very pronounced 
in the case of the fcc(100) structure. This holds not only for thin films where some 
oscillations are present due to the finite film thickness, but also extends to the limit 
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Figure 3. Surface-, center-, and mean magnetization (mi, m c , m) as a function of 
the film thickness d for three different temperatures T = 0.1 Tc, T = 0.8 Tq, and 
T = 0.9 T c ; (a) fcc(100); (b) fcc(lll). Further parameters: U = 50eV, t = -0.25 eV, 
n = 1.6. 



d — > oo where the two surfaces are well separated and do not interact. The oscillations 
as a function of d which are present for the fee (100) structure get damped for higher 
temperatures. One can see from Fig. [3] that the center layer magnetization m c for thick 
films (d > 10) is in good agreement with the corresponding fee bulk calculation [49 . 

The magnetization profile for both film geometries is plotted for d = 10 in Fig. ||. 
Here again the magnetizations of the fcc(100) film show a pronounced layer dependence 
while they are very close to the bulk value from the second layer on in the case of the 
fcc(lll) film geometry. The magnetization profiles are similar to the ones obtained in 
PH for Cu/Ni/Cu sandwiches calculated within a single-site spin-fluctuation theory. 
However, within the present approach the deviation from the bulk magnetization is 
enhanced close to the Curie temperature for the fee (100) structure (Fig. |). Note 
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Figure 4. Magnetization profile of a ten layer film for three different temperatures 
T = 0.1 T c (solid line), T = 0.8T C (dotted line), and T = 0.9 T c (dashed line). Further 
parameters: U — 50 eV, t = — 0.25eV, n = 1.6. 



that a similar trend to a reduced surface magnetization is also found within localized 
spin models. However, for the uniform Heisenberg model without a layer-dependent 
anisotropy contribution, the layer magnetizations necessarily increase monotonously 
from the surface to the central layer P2, PJ. 

To understand the magnetic behavior on a microscopic basis we will, in the 
following, discuss the temperature-dependent electronic structure of the thin film 
systems. For a five layer fee (100) film the spin- and layer-dependent spectral density 
at the gamma f point and the quasiparticle density of states are plotted in Fig. |. 
There appear two correlation induced band-splittings in the quasiparticle spectrum: 
Due to the strong Coulomb interaction the spectrum splits into a low and a high 
energy subband ( "Hubbard bands" ) which are separated by an energy of the order U. 
Besides this so-called "Hubbard splitting" that is present for all temperatures there is 
an additional exchange splitting in majority (p =f ) and minority (a =|) spin direction 
for temperatures below Tc- In the lower subband the electron mainly hops over empty 
sites, while in the upper subband it hops over lattice sites that are already occupied by 
another electron with opposite spin. The corresponding weights of the subbands scale 
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with the probability of the realization of these two situations. In the strong coupling 
limit the scaling is given by (1 — n a _ a ) and n Q ,_ cr for the lower and upper subband, 
respectively. Since the total band occupation (n = 1.6) considered here is above half 
filling, the chemical potential // is located in the upper subband while the lower subband 
is completely filled. 

Starting from the Curie temperature the evolution of the quasiparticle spectrum 
with decreasing temperature is dominated by two distinct correlation effects. Both are 
driven by an increasing spin- asymmetry in the bandshift B a _ a . Firstly the centers of 
gravity of the majority and minority subbands move apart with decreasing temperature 
(Stoner-type behavior). Secondly there is a strong spin-dependent transfer of spectral 
weight between the lower and the upper subbands according to the above mentioned 
scaling, which results in spin- and temperature-dependent widths of the respective 
subbands. This behavior can also be seen in detail in Fig. || and Fig. [7| where the 
quasiparticle bandstructure of the surface and central layer is plotted for a five layer 
fcc(lOO) and fcc(lll) film, respectively. Here, only the upper subbands are shown. 
While the centers of gravity of the upper a =[ subbands are shifted to higher energies 
for decreasing temperatures, the lowest excitation peak in the spectral density at T is 
even lowered due to the increasing bandwidth. On the other hand the width of the 
upper a =] subband decreases. The interplay of these two correlation effects leads to 
an inverse exchange splitting at the lower edge of the upper subband near to the T 
point. The corresponding quasiparticle density of states is, however, very small. Note, 
that for the same reason the position of the central peak of the upper subband is almost 
spin- and temperature- independent. This behavior holds for both film structures for 
k- vectors not too far away from f. 

With help of the quasiparticle bandstructure given in Fig. ^] and Fig. [7| the exchange 
splitting between majority and minority spin direction can be analyzed in more detail. 
For both film structures the exchange splitting is wavevector-dependent. It is strongest 
near M for fcc(lOO) and between M and K for fcc(lll). Contrary to the fcc(lOO) 
structure all layers are fully polarized at T = OK for the fcc(lll) film. In the case 
of ferromagnetic saturation the exchange splitting can be estimated [48] to be at 
most (1 — n a {)[At — B a ±], where (1 — n a i)B a i is the effective bandshift between the 
centers of gravity of the upper quasiparticle subbands. For strong Coulomb interaction 
n a i(l—n a i)B a i is proportional to the kinetic energy J|7| of the a = j electrons in the a-th 
layer. Note that the kinetic energy of the a =\ electrons vanishes for the ferromagnetic 
saturated state since the a =\ band is completely filled. We want to point out that 
these results strongly contrast the findings of Hartree-Fock theory where the exchange 
splitting is wavevector-independent and proportional to m a U leading to substantially 
higher Curie temperatures compared to the MAA. The temperature-dependence of the 
electronic structure within the MAA is completely different to the Stoner picture of 
ferromagnetism. 

Let us discuss the quasiparticle lifetime which corresponds to the width of the 
quasiparticle peaks. From the spectral density (Figs. ||, ||, |7|) one can clearly read off 
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that the lifetime of the quasiparticles is strongly spin- and temperature-dependent. For 
low temperatures the upper minority spectrum is sharply peaked which indicates long 
living quasiparticles. This is due to the fact that in the ferromagnetic saturated state a 
a =1 electron does meet a a =] electron at any lattice site and thus effectively does not 
perform any scattering process. The width of the a —\ quasiparticle peaks, however, is 
broadened for decreasing temperature. Thus in the majority spectrum the quasiparticle 
lifetime decreases for increasing magnetization. What concerns the lower subbands (see 
Fig. H) the respective spectrum is strongly damped and the different excitations due to 
the five layer structure are almost indistinguishable. 

For given spin and wavevector the positions of the quasiparticle peaks are layer- 
independent. In principle, their number corresponds to the number of layers of the film. 
However, due to symmetry some peaks are left out for certain layers. For thicker films 
the different peaks move closer together as their number increases until they build a 
continuum for d — > oo which corresponds to the projection of the three dimensional 
bandstructure onto the surface Brillouin zone. Between M and X for fcc(lOO) and at 
K for fcc(lll) the different peaks merge together due to vanishing interlayer hopping 
( 7 j_(k) = 0). ' 

In Figs. ||, |7| the QDOS of the surface and central layer are shown additionally. 
The Van Hove singularities resulting from the different branches of the quasiparticle 
dispersion are clearly visible. There are sharp Van Hove singularities in the minority 
spectrum while they are broadened for the majority spin direction because of the finite 
widths of the a =] quasiparticle peaks due to the quasiparticle damping. 

Finally we want to stress that the results presented above do not depend on the 
size of the Coulomb interaction U as long as U is chosen from the strong-coupling region 
(U ^> W). Contrary to Hartree-Fock theory, all magnetic key quantities like the Curie 
temperature and the exchange splitting saturate as a function of U. On the other hand, 
although the MAA was optimized with respect to the strong coupling limit we believe 
that, at least qualitatively, the correlation effects in the spin-, layer-, and temperature- 
dependent electronic structure are valid down to intermediate Coulomb interaction as 
well. 

6. Conclusion 

For the investigation of spontaneous ferromagnetism and electron correlation effects 
in thin itinerant-electron films we have applied a generalization of the modified alloy 
analogy (MAA) to the single-band Hubbard model with reduced translational symmetry. 
The MAA is based on the alloy analogy concept and is optimized with respect to 
correct strong coupling behavior Jill |45j. Within the MAA the actual type of the 
underlying alloy is not predetermined but has to be determined selfconsistently. In 
this sense the MAA is able to account for the itineracy of the —a electrons which 
are considered as strictly "frozen" in the conventional alloy analogy (AA). In the 
paramagnetic phase MAA and AA are almost identical. However, contrary to the AA 
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Figure 7. Spin and layer dependent spectral density and quasiparticle density of states 
of a five layer fcc(lll) film for three different temperatures T = 0.17c, T = 0.9 Tc, 
and T = Tc- Only the upper Hubbard band is shown, a = 1: surface layer; a = 3: 
central layer. The chemical potential is located at zero energy. Further parameters: 
U = 50 cV, t = -0.25 eV, n = 1.6. 
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spontaneous ferromagnetic order is possible for special parameter constellations within 
the MAA. With help of the MAA the interplay of magnetism and quasiparticle damping 
effects can be studied in a natural way. 

For an fcc(100) and an fcc(lll) film geometry the layer-dependent magnetizations 
have been discussed as a function of temperature as well as film thickness. The 
magnetization in the surface layer is found to be reduced with respect to the inner 
layers for all thicknesses and temperatures considered. While this reduction is weak for 
fcc(lll) films it is pronounced in the case of an fcc(100) geometry. The effect of the 
surface is considerably stronger for fee (100) films due to the higher percent of missing 
nearest neighbor atoms. The reduction of the surface layer magnetization is not to be 
expected within an Hartree-Fock type approach (Stoner criterion) to the Hubbard film 
being, therefore, a genuine effect induced by strong electron correlations. 

The magnetic behavior of the thin film systems can be microscopically understood 
by means of the spin- layer- and temperature- dependent quasiparticle bandstructure 
and the corresponding quasiparticle density of states. There appear two correlation 
induced band splittings in the quasiparticle spectrum. Besides the Hubbard 
splitting there is an additional exchange splitting for temperatures below Tq- The 
demagnetization process as a function of temperature is dominated by two distinct 
correlation effects: A Stoner-like shift in the centers of gravity of the majority and 
minority subbands together with a strong spin-dependent transfer of spectral weight 
between the upper and lower subbands. An interplay of these two effects results in Curie 
temperatures far below the corresponding Hartree-Fock values. The exchange splitting is 
found to be strongly wavevector-dependent and is substantially different for the various 
quasiparticle branches in the bandstructure. The widths of the quasiparticle peaks 
that correspond to the quasiparticle lifetime exhibit a strong spin- and temperature- 
dependence. For T = K the minority-spin quasiparticle peaks are sharply peaked 
while the majority-spin spectrum is substantially broadened. 

Clearly the degeneracy of the 3d-bands has to be included if a direct comparison 
to the experiment is intended. Within the present scheme this could be achieved by a 
similar approach as presented in |33] which is planed for the future. However we believe 
the correlation effects found here to be important within a generalized Hubbard model 
as well. In this work we have exclusively focused on purely ferromagnetic films. In 
addition one can examine within the same theory a phase with antiferromagnetic order 
between the layers. We expect such a situation to exist close to half-filling (n = 1) 
and for intermediate values of the Coulomb interaction. Further the influence of a 
non-magnetic top layer on the magnetic behavior of thin films can be investigated. 
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